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Abstract 


Background: Vaccine mandates and vaccine passports (VMVP) for SARS-CoV-2 are thought to 
be a path out of the pandemic by increasing vaccination through coercion and excluding 
unvaccinated people from different settings because they are viewed as being at significant risk 
of transmitting SARS-CoV-2. While variants and waning efficacy are relevant, SARS-CoV-2 
vaccines reduce the risk of infection, transmission, and severe illness/hospitalization in adults. 
Thus, higher vaccination levels are beneficial by reducing healthcare system pressures and 
societal fear. However, the benefits of excluding unvaccinated people are unknown. 


Methods: A method to evaluate the benefits of excluding unvaccinated people to reduce 
transmissions is described, called the number needed to exclude (NNE). The NNE is analogous 
to the number needed to treat (NNT=1/ARR), except the absolute risk reduction (ARR) is the 
baseline transmission risk in the population for a setting (e.g., healthcare). The rationale for the 
NNE is that exclusion removes all unvaccinated people from a setting, such that the ARR is the 
baseline transmission risk for that type of setting, which depends on the secondary attack rate 
(SAR) typically observed in that type of setting and the baseline infection risk in the population. 
The NNE is the number of unvaccinated people who need to be excluded from a setting to 
prevent one transmission event from unvaccinated people in that type of setting. The NNE 
accounts for the transmissibility of the currently dominant Delta (B.1.617.2) variant to estimate 
the minimum NNE in six types of settings: households, social gatherings, casual close contacts, 
work/study places, healthcare, and travel/transportation. The NNE can account for future 
potentially dominant variants (e.g., Omicron, B.1.1.529). To assist societies and policymakers in 
their decision-making about VMVP, the NNEs were calculated using the current (mid-to-end 
November 2021) baseline infection risk in many countries. 


Findings: The NNEs suggest that at least 1,000 unvaccinated people likely need to be excluded 
to prevent one SARS-CoV-2 transmission event in most types of settings for many jurisdictions, 
notably Australia, California, Canada, China, France, Israel, and others. The NNEs of almost 
every jurisdiction examined are well within the range of the NNTs of acetylsalicylic acid (ASA) 
in primary prevention of cardiovascular disease (CVD) (= 250 to 333). This is important since 
ASA is not recommended for primary prevention of CVD because the harms outweigh the 
benefits. Similarly, the harms of exclusion may outweigh the benefits. These findings depend on 
the accuracy of the model assumptions and the baseline infection risk estimates. 


Conclusions: Vaccines are e beneficial, but the eee NNEs suggest that excluding unvaccinated 
ogli : nsmissions in many E across the eee 





This i 15  becausd | | 
e EE most types of settings s since current baseline transmission 
risks are negligible. Consideration of the harms of exclusion is urgently needed, including 
staffing shortages from losing unvaccinated healthcare workers, unemployment/unemployability, 
financial hardship for unvaccinated people, and the creation of a class of citizens who are not 
allowed to fully participate in many areas of society. 
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INTRODUCTION 

A key global priority for the pandemic caused by severe acute respiratory syndrome 
coronavirus 2 (SARS-CoV-2) is to return to normality as quickly as possible. This is because 
many societies are losing their political, economic, psychological, and social tolerance for 
lockdowns and restrictions, especially since the development of SARS-CoV-2 vaccines. There is 
growing pressure to implement vaccine mandates and vaccine passports (VMVP) domestically 
and for international travel as a “ticket” for full participation in post-pandemic life. 
Consequently, more governments, healthcare organizations, educational institutions, and 
businesses are adopting VMVP. VMVP are thought to be a path out of the pandemic by (i) 
increasing vaccination levels through coercion and (ii) excluding those who remain unvaccinated 
from many types of settings (e.g., work, social gatherings). A central legal justification for 
exclusion is that unvaccinated people are viewed as being at significant risk of transmitting 
SARS-CoV-2, which is potentially the most serious externalized risk of the choice to remain 
unvaccinated.' 

The benefits of vaccination in response to VMVP are clear. Meta-analyses show that 
SARS-CoV-2 vaccines are effective at reducing the risk of infection, transmission, and severe 
illness/hospitalization in adults.?3 Therefore, while variants and waning efficacy may moderate 
these benefits, higher vaccination levels are beneficial because they will reduce the burden on 
healthcare systems and reduce societal fear. On the other hand, the benefits of excluding 
unvaccinated people to reduce SARS-CoV-2 transmissions has not been evaluated quantitatively. 
This evaluation is a critical first step to determining the costs vs. benefits of VMVP. It is also 
urgently needed because an unknown percentage of people will not get vaccinated despite 


VMVP. This means those unvaccinated people will be excluded from many areas of modern life, 
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such as many places of work, leisure, transportation, and education. As of December 2021, this is 
already happening in many jurisdictions globally. 

Therefore, we performed a pre-planned analysis of a larger study of VMVP 
(CRD42021292263). The purpose of this analysis was to describe a method to evaluate the 
benefits of excluding unvaccinated people to reduce transmissions analogous to the number 
needed to treat (NNT),* which we call the number needed to exclude (NNE). The NNT is 1 
divided by the absolute risk reduction (ARR). The NNT is the number of people who need to 
receive a treatment to prevent one outcome (e.g., myocardial infarction). The difference is that 
the ARR of the NNE is the baseline transmission risk in the population for a given type of setting 
(NNE=1/ARR). The rationale for the NNE is that exclusion removes all unvaccinated people 
from a setting (e.g., healthcare), such that the ARR is the baseline transmission risk in that type 
of setting. The ‘treatment’ is the public health intervention of exclusion via VMVP. 

The baseline transmission risk is the probability in the general population that an 
unvaccinated person is currently infected and transmits SARS-CoV-2 in a given type of setting. 
In other words, it is the current probability in the general population of a transmission event in a 
setting from an unvaccinated person. The reciprocal of this probability is the NNE. The baseline 
transmission risk is estimated by taking the combined probability of the secondary attack rate 
(SAR) typically observed in that type of setting and the baseline infection risk in the general 
population. The combined probability is needed to estimate this probability because an 
unvaccinated person must be infected first before they can transmit SARS-CoV-2. Therefore, the 
NNE is the number of unvaccinated people who need to be excluded from a setting to prevent 
one transmission event from unvaccinated people in that type of setting. The NNE pertains to 


one transmission event, which may include one or more secondary infections, because the SAR 
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is the infection risk among the contacts of the unvaccinated index case. The NNE concerns the 
direct transmission caused by the unvaccinated person, not subsequent generations of 
transmission. 

The baseline infection risk is the current point-prevalence of infectious cases in the 
general population. The numerator is the total number of new and existing infectious cases, and 
the denominator is the population size. It is the ‘baseline’ risk because it is the current risk of 
getting infected with SARS-CoV-2 in a population. It is the estimated risk that an unvaccinated 
person gets infected. Point-prevalence is a more appropriate metric of baseline infection risk than 
incidence or period prevalence for three reasons. First, the population-level risk of getting 
infected depends not just on new infectious cases, but existing ones too. Second, while risk over 
time is important, people are primarily concerned about the risk of getting infected with SARS- 
CoV-2 right now, not, for example, over the last 3 months. Second, incidence depends on time at 
risk. Shorter periods will have a lower incidence than longer periods. There is no non-arbitrary 
way of defining the “correct” time at risk. Point-prevalence does not have this arbitrariness 
because it uses the current prevalence. Like the NNT, time is implicit in the NNE since it relates 
to the period over which the risk is measured. The NNE relates to the current risk because it is 
based on the current point-prevalence of infectious cases and thus the current baseline 
transmission risk. In practice, the risk would be over one day given how point-prevalence data is 
typically reported. Time is also incorporated into the NNE when one calculates the NNEs over 
time using changing baseline infection risks. 

Critically, societies need to be mindful of the baseline infection risk because when it is 
low, baseline transmission risks are also low, which exponentially increases the NNE (Figure 1). 


Figure 1 shows that the transmission risk reduction gained from excluding unvaccinated people 
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becomes negligible when baseline transmission risks are low, especially below 0.40% to 0.30%. 
This is the range of the ARRs of acetylsalicylic acid (ASA) in primary prevention of 
cardiovascular disease (CVD) (NNTs > 250 to 333)?” Notably, ASA is not recommended for 
primary prevention in all adults because the negligible ARRs do not outweigh the risk of harm.’ 


Similarly, the harms of exclusion via VMVP may outweigh the benefits when baseline infection 


risks are low. 
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Figure 1. Relationship between baseline transmission risk and the number needed to exclude (NNE). The 
number needed to treat (NNT) has the same scale as the NNE (NNT=1/ARR). The NNTs of other 
interventions are shown for comparison.**°*!! The ARRs of the influenza and SARS-CoV-2 vaccines are 
for reducing the risk of infection. The ARR of the HPV vaccine is for preventing any cervical pre-cancer. 
ARR=absolute risk reduction. MI=myocardial infarction. HPV=human papillomavirus. 
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The NNE was calculated for six types of settings: (1) households, (ii) social gatherings 
(e.g., meals, conversations, social gatherings with friends/family), Gi) casual close contacts 
(e.g., public areas or buildings), (iv) work or study places, (v) healthcare, and (vi) 
travel/transportation. The NNEs accounted for the transmissibility of the Delta (B.1.617.2) 


variant because it currently predominates in most countries. The NNE can account for the 
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increased transmissibility of future potentially dominant variants (e.g., Omicron, B.1.1.529). 
First, the NNEs for each type of setting were modeled as a function of simulated baseline 
infection risks (0.10% to 100%). Second, to assist societies and policymakers in their decision- 
making about excluding unvaccinated people via VMVP, the NNEs were calculated for each 
type of setting using the current (mid-to-end November 2021) baseline infection risk in many 
countries around the globe. 
METHODS 

The NNE for the ith type of setting was modeled as follows: 


1 1 
` ARR; (SARy, X CFpeuta) X BIR 


NNE; 
The ARR of excluding unvaccinated people from the ith setting is the baseline transmission risk 
for that type of setting. The baseline transmission risk is estimated by multiplying robust 
estimates of the wild-type SAR for the ith setting (SAR) by a correction factor for the 
increased transmissibility of the Delta variant over the wild-type (CFpej¢q) and the baseline 
infection risk (BIR) in the population. CFpelta is necessary because our robust SARs were 
derived from meta-analyses of studies from 2020 when the wild-type form of SARS-CoV-2 
predominated, whereas the Delta variant currently predominates in most counties. We needed to 
use the wild-type SARs because studies of the SARs of variants occurred during the vaccine 
period (2021), which would add a confounding factor to the model. CFperta was 1.97 in the 
model because a global analysis of effective reproduction numbers found that the Delta variant 
was 97% more transmissible than the wild-type.'? The model can account for future potentially 
dominant variants (e.g., Omicron, B.1.1.529) as data emerges regarding its transmissibility over 
the wild-type by changing the value of CF. The corrected robust SARs (SARw, X CFpeita) must 
be multiplied by the BIR because SARs are the transmission risks among infected people, not the 
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general population. Estimating the transmission risk in the general population requires 
calculating the combined probability of infection and transmission because a person must be 
infected first before they can transmit SARS-CoV-2. To model how the NNE for the ith setting 
changes as a function of the baseline infection risk, we simulated values of BIR ranging from 
0.10% to 100% and plotted the results. 

Robust estimates of the wild-type SARs (SARy,) were obtained by calculating the mean 
and 95% confidence intervals (CIs) of meta-analytic estimates of the SARs for the six types of 
settings mentioned above. Social gatherings are high intensity contacts (e.g., gatherings of 
friends/family), whereas casual close contacts are lower intensity contacts (e.g., contact in a 
public area or building). The 95% CIs of the NNEs were calculated from the upper and lower 
limits of these robust SARs. To identify these meta-analyses, we searched PubMed from 
inception to November 26, 2021, with no language restrictions. We used the search terms 
"SARS-CoV-2" or "COVID-19" combined with “secondary attack*” and “meta-analysis”. Two 
authors (AP, BH) screened the titles/abstracts to determine potential eligibility. Full-texts of 
potentially eligible reports were retrieved and assessed to determine eligibility. Data were 
independently extracted. Disagreements at either the screening, full-text review, or data 
extraction stages were resolved through consensus between the two authors. Reports were 
included if they were a meta-analysis of pre-vaccine studies (2020) of SARs and estimated the 
SAR for at least one of the six types of settings. Meta-analyses of studies during the vaccine 
period (2021) were excluded to remove the confounding influence of vaccination on SARs. 
Random effect model data were extracted whenever possible. 13 citations were identified from 


the search. 1 meta-analysis! was identified in our pilot study.* The meta-analytic estimates of 
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the SARs were extracted from the 7 meta-analyses meeting the inclusion criteria.!3-1? All 
analyses were conducted in R.”° 

The NNE model makes four assumptions. First, we assume the Delta variant correction 
factor (CFpeitq=1.97) is accurate. We checked this assumption by comparing the model 
predictions with the observed SARs of the Delta variant in other studies (see Results). Second, 
the model assumes the SARs for each setting type are relatively consistent because of systematic 
differences in the intensity of contact (proximity, duration of contact, use of precautions, e.g., 
masking). Lower contact intensity settings should have lower SARs than higher contact intensity 
settings. For instance, we assume that households will consistently have higher SARs than 
healthcare settings because the intensity of contact is higher, which is suggested by the fact that 
households are a major driver of SARS-CoV-2 transmissions.'© We supported this assumption by 
utilizing meta-analytic estimates in order to arrive at robust estimates of the central tendency and 
distribution of SARs for each type of setting. 

Third, the model assumes the baseline infection risk is relatively stable across individuals 
in a population. Given that there will inevitably be local variation in risk, we addressed this 
assumption by modeling a range of baseline infection risks from 0.10% to 100% and plotted the 
95% CIs around the NNEs based on the 95% CIs of the corrected robust SARs. Similarly, we 
calculated the 95% CIs around the NNEs for each jurisdiction we examined (see below) using 
the 95% CI of the baseline infection risk in the jurisdiction. 

Fourth, the model assumes (1) unvaccinated people have no natural immunity and (ii) the 
contacts of infected unvaccinated index cases have no vaccine or natural immunity. This is 
because the SARs were derived from the pre-vaccine period. It was not possible to know from 


the meta-analyses to what extent the studies included participants with a history of prior infection 
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with SARS-CoV-2. However, it is reasonable to assume that the vast majority had no natural 
immunity in these studies. This is because a global meta-analysis of studies from 2020 estimated 
the median seroprevalence of SARS-CoV-2 antibodies in the general population across countries 
was 4.5% (IQR: 2.4%, 8.4%).21 This is important because populations with higher levels of 
vaccine and/or natural immunity will have lower SARs and thus lower baseline transmission 
risks, regardless of the baseline infection risk in the population. Accordingly, higher levels of 
vaccine/natural immunity in a population will increase the NNE, independent of the current 
baseline infection risk. By assuming this, the corrected robust SARs (SARw, X CFpeita) are the 
worst-case scenario because they are higher than what would be predicted in a population with 
some level of vaccine/natural immunity. This is a strong conservative assumption since many 
jurisdictions already have high levels of vaccine/natural immunity (as of December 2021). 
Vaccine/natural immunity in a population will reduce the denominator of the NNE because they 


will reduce the value of SARw,, which will reduce the baseline transmission risk, and thus 


increase the NNE. Therefore, the calculated NNEs likely represent the minimum number of 
unvaccinated people needed to exclude to prevent one SARS-CoV-2 transmission event in each 
type of setting. 

Research supports this impact of vaccine/natural immunity on the baseline transmission 
risk in a population. A large Belgian study of the period of Alpha (B.1.1.7) predominance found 
that natural immunity had the same effectiveness (~50-60%) for reducing the SAR among 
unvaccinated close contacts as full vaccination with BNT162b2 and mRNA-1273.” It was also 
the same (~75-85%) for contacts with natural immunity vs. vaccine immunity who were exposed 
to unvaccinated index cases. A large Dutch study during the same period found that the relative 


risk reduction (RRR) is 65% when unvaccinated index cases had contact with fully vaccinated 
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(SAR=11%) vs. unvaccinated (SAR=31%) household members.” In summary, this research 
shows that if an unvaccinated index case has natural immunity or their contacts have either 
vaccine or natural immunity, the SAR is reduced for the contacts of infected unvaccinated 
people. Therefore, vaccine/natural immunity in a population will reduce the baseline 
transmission risks in that population for any given setting, and thus increase the NNE. 

One can incorporate the effect of vaccine/natural immunity into the NNE model by 
adding additional correction factors to the corrected robust SARs (SAR, X CFpetta). These 
correction factors can be estimated using the RRRs of vaccine/natural immunity on SARs in 
unvaccinated index cases and their contacts, each weighted by the proportion of the population 
fully vaccinated and with natural immunity. However, we opted not to do this for two reasons. 
First, we wanted to make as few assumptions as possible so that the model remained close to the 
data. Second, we wanted to produce conservative estimates of the NNEs for each setting. 

The NNEs were also calculated using estimates of current baseline infection risks and 
their 95% Cls for a sample of countries across the globe. The same computations above were 
used, except the 95% Cls of the NNEs were calculated from the upper and lower limits of the 
baseline infection risks. Estimates of the current point-prevalence of infectious cases in these 
jurisdictions were taken from a global database (https://decision-support-tools.com/) developed 
and regularly updated by Defence Research and Development Canada (DR DO), an agency of the 
Canadian Department of National Defence. Its estimation methods are detailed online. Australia 
had data on its six states, the Australian Capital Territory, and Northern Territory. China had data 
on its provinces, municipalities, and autonomous regions except Tibet. For both countries, these 
data were extracted and means were used to estimate national-level risks. No data for Sweden 


was available for the time point we extracted. Regional data for Canada and the United States of 
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America (USA) were available, extracted, and means were used to estimate state/provincial-level 
and national-level risks. We extracted the most recent dataset posted online on November 261, 
2021. 71% of the dataset was from November 21 to 25", 2021, 29% were from November 7" to 
11", 2021. Thus, point-prevalence data were from mid-to-end November 2021. 
RESULTS 

Table | displays the wild-type vs. corrected robust SARs for each type of setting. The 
predicted household SAR for the Delta variant (32.59%, 95% CI: 27.04%, 38.14%) was 26% 
higher than the observed SARs pooled across recent studies of unvaccinated index cases infected 
with the Delta variant who had contact with primarily household members (25.90%, 95% CI: 
19.02%, 32.78%).°*° In other words, the Delta variant correction factor (1.97) based on the 
effective reproduction numbers was an accurate but conservative assumption because it 
increased the SARs more than what has been observed in the real-world. This additional 


conservative assumption of the model is another reason why the NNEs are likely the minimum 








NNEs. 
Table 1. Robust estimates of the secondary attack rates to model the NNE for each setting 
Corrected robust SARs 
Robust SARs (Delta variant: 
Type of setting (Wild-type: SARw,) (95% CD SARw, X CF peita) 95% CD 
Households 16.54% (13.73%, 19.36%) 32.59% (27.04%, 38.14%) 
Social gatherings 5.93% (5.87%, 6.00%) 11.69% (11.56%, 11.82%) 
Casual close contacts 1.55% (0.86%, 2.24%) 3.05% (1.70%, 4.40%) 
Work/study places 1.47% (0.96%, 1.98%) 2.89% (1.88%, 3.89%) 
Healthcare 1.50% (0.12%, 2.88%) 2.96% (0.24%, 5.67%) 
Travel/transportation 2.23% (0.00%, 4.95%) 4.40% (0.00%, 9.75%) 





95% Cl=95% confidence interval. Brackets are the 95% confidence intervals (lower limit, 
upper limit). SAR=secondary attack rate. NNE=number needed to exclude. Note, the CIs were 
widest for travel/transportation and because the lower limit crossed zero it was set to zero. 

CF petta=1.97. See Methods for details. 


Figure 2 plots the NNEs for the six types of settings using simulated baseline infection 
risks. These NNEs can be interpreted in terms of the NNTs of ASA in primary prevention of 


CVD (Figure 1). NNEs within the range of the NNTs of ASA in primary prevention of CVD (= 
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250 to 333) are considered ‘high’ because it means the benefits of excluding unvaccinated people 
to reduce transmission risks are negligible, much like the ARRs of ASA are negligible in 
preventing myocardial infarction and strokes in people with no history of CVD.?” These 
negligible ARRs result in high NNTs and are the reason why ASA is not recommended for 
primary prevention of CVD in all adults since the benefits do not outweigh the risk of harm.’ As 
seen in Figure 2, for household settings, NNEs start to become high when baseline infection risks 
are < 1%. For social gatherings, NNEs start to become high when baseline infection risks are < 
3%. For casual close contacts, work/study places, and healthcare settings, NNEs start to become 
high when baseline infection risks are < 10%. For travel/transportation settings, NNEs start to 
become high when baseline infection risks are < 5-10%. Figure 2 shows that when baseline 
infection risks are below 1%, baseline transmission risks become very small, such that thousands 
of unvaccinated people likely need to be excluded from most types of settings to prevent one 
SARS-CoV-2 transmission event. High NNEs mean that, when the current point-prevalence of 
infectious cases in the population is below these thresholds, unvaccinated people are likely not at 


significant risk — in absolute terms — of transmitting SARS-CoV-2 to others. 
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Figure 2. Modeling the number of unvaccinated people needed to exclude (NNE) from different types of settings to prevent one SARS-CoV-2 
transmission event in that setting. The NNEs are plotted as a function of simulated baseline infection risks (0.10% to 100%) and corrected for the 
increased transmissibility of the Delta variant. To better visualize the NNEs at different baseline infection risks, each setting is plotted from 0.10% 
to 0.50%, 1% to 5%, and 0% to 100%. 


medRxiv preprint doi: https://doi.org/10.1101/2021.12.08.21267162; this version posted December 11, 2021. The copyright holder for this 
preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. 


It is made available under a CC-BY-NC-ND 4.0 International license . 


Table 2 displays the NNEs for many jurisdictions from the DRDC database estimates of 
the baseline infection risk as of mid-to-end November 2021. Many of these jurisdictions have 
implemented some form of VMVP resulting in the exclusion of unvaccinated people from 
different areas of these societies. NNEs are color-coded by the lower limit of the 95% CI: yellow 
(= 1,000 and < 10,000), orange (= 10,000 and < 100,000), red (= 100,000 and < 1,00,000), 
and black (> 1,000,000). NNEs > 1,000,000 were truncated to > 1M. The NNEs of almost 
every jurisdiction and setting type are well within the range of the NNTs of ASA in primary 
prevention of CVD (= 250 to 333). The NNEs suggest that at least 1,000 unvaccinated people 
likely need to be excluded to prevent one SARS-CoV-2 transmission event in most types of 
settings for many jurisdictions, notably Australia, California, Canada, China, France, Israel, and 


others. Given the conservative assumptions of the model, these are likely the minimum NNEs. 
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Table 2. Global estimates of the number of unvaccinated people needed to exclude (NNE) to prevent one SARS-CoV-2 transmission event 





Estimated baseline 
infection risks as of mid- 


Estimated number of unvaccinated people needed to exclude (NNE) 
to prevent one SARS-CoV-2 transmission event in each type of setting (95% CI) 




































to-end November 2021+ Travel/ 
Jurisdiction (95% CD Households Social gatherings Casual close contacts ` Work/study places Healthcare transportation 
Argentina 0.0798% 3,843 10,716 41,019 43,349 42,386 28,468 
(0.0663%, 0.1091%) (2,813, 4,626) (7,842, 12,898) (30,019, 49,373) (31,724, 52,178) (31,019, 51,019) (20,834, 34,266) 
Australia 0.2391% 1,283 3,578 13,695 14,473 14,151 9,504 
(0.135%, 0.3425%) (896, 2,273) (2,498, 6,338) (9,562, 24,262) (10,106, 25,640) (9,881, 25,071) (6,637, 16,839) 
Australian Capital 0.2908% 1,055 2,942 11,262 11,901 11,637 7,816 
Territory (0.1495 %, 0.5032%) (610, 2,052) (1,700, 5,722) (6,508, 21,904) (6,878, 23,149) (6,725, 22,634) (4,517, 15,202) 
New South Wales 0.0903 % 3997 9,471 36,255 38,315 37,464 25,162 
(0.0614%, 0.1252%) (2,450, 4,998) (6,832, 13,934) (26,153, 53,341) (27,639, 56,371) (27,025, 55,119) (18,151, 37,020) 
Northern Territory 0.3234% 949 2,645 10,125 10,700 10,463 7,027 
(0.1742%, 0.6011%) (510, 1,761) (1,423, 4,910) (5,448, 18,795) (5,757, 19,863) (5,629, 19,421) (3,781, 13,044) 
Queensland 0.0003% > IM > IM > IM > IM > IM > IM 
(0.0001 %, 0.001%) (320,491, = 1M) (893,569, > 1M) (= 1M, = 1M) (> 1M, = IM) (2 1M, = IM) (= 1M, = IM) 
South Australia 0.0002% > IM > IM > IM > IM = IM > IM 
(0.0001 %, 0.0005%) (606,215, = 1M) (= 1M, = 1M) (= 1M, = 1M) (> 1M, = IM) (> 1M, = IM) (= 1M, = 1M) 
Tasmania 0.0009% 337,842 941,943 > IM > IM > IM > IM 
(0.0003%, 0.002%) (157,044, 902,176) (437,857, = 1M) (= 1M, = 1M) (= 1M, = 1M) (= 1M, = 1M) (= 1M, = 1M) 
Victoria 1.2068% 254 709 2,714 2,868 2,804 1,883 
(0.694%, 1.5058%) (204, 442) (568, 1,233) (2,175, 4,719) (2,298, 4,987) (2,247, 4,876) (1,509, 3,275) 
Western Australia 0.0004% 817,673 > IM > IM > IM > IM > IM 
(0.0002 %, 0.001%) (292,777, = 1M) (816,298, = IM) (= 1M, = 1M) (> 1M, = IM) (= 1M, = 1M) (= 1M, = 1M) 
Brazil 0.458% 670 1,868 7,151 TSS 7,389 4,963 
(0.302%, 0.7268%) (422, 1,016) (1,177, 2,833) (4,506, 10,846) (4,762, 11,462) (4,656, 11,207) (177527) 
Canada 1.1221% 273 762 2,919 3,085 3,016 2,026 
(0.7015%, 1.7193%) (178, 437) (498, 1,220) (1,905, 4,669) (2,013, 4,934) (1,968, 4,824) (1,322, 3,240) 
Alberta 0.6645% 462 1,288 4,929 5,209 5,093 3,421 
(0.356%, 1.2135%) (253, 862) (705, 2,403) (2,699, 9,200) (2.552 0 22) (2,789, 9,506) (1,873, 6,385) 
British Columbia 1.0078% 280 779 2,983 3,153 3,083 2,070 
(0.7027%, 1.5806%) (194, 437) (541, 1,217) (2,072, 4,660) (2,190, 4,925) (2,141, 4,816) (1,438, 3,235) 
Manitoba 2.2159% 138 386 1,478 1,562 (Eer 1,026 
(1.3999%, 3.2252%) (95, 219) (265, 611) (1,015, 2,339) (1,073, 2,472) (1,049, 2,417) (705, 1,624) 
New Brunswick 1.2565% 244 681 2,606 2 755 2,693 1,809 
(0.8438%, 1.7817%) (172, 364) (480, 1,014) (1,838, 3,881) (1,943, 4,102) (1,899, 4,011) (1,276, 2,694) 
Newfoundland and 0.034% 9,030 25,178 96,381 101,857 99,593 66,891 
Labrador (0.0147%, 0.0639%) (4,804, 20,878) (13,395, 58,212) (51,275, 222,832) (54,188, 235,492) (52,984, 230,259) (35,586, 154,652) 
Northwest Territories 2.0711% 148 413 1,581 1,671 1,634 1,097 
(1.0044%, 3.4395%) (89, 305) (249, 852) (952, 3,260) (1,006, 3,446) (984, 3,369) (661, 2,263) 
Nova Scotia 0.3008% 1,020 2,844 10,888 11,507 11,251 7,557 
(0.1791%, 0.4943%) (621, 1,714) (1,731, 4,778) (6,625, 18,290) (7,001, 19,329) (6,846, 18,899) (4,598, 12,694) 
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Nunavut 0.0000% > IM 
(0.0000%, 0.0011%) (269,662, = 1M) 
Ontario 0.9107% 337 
(0.5362%, 1.4796%) (207, 572) 
Prince Edward Island 0.7575 % 405 
(0.4458%, 1.0521%) (292, 688) 
Quebec 1.611% 190 
(1.0858%, 2.3744%) (129, 283) 
Saskatchewan 0.6956% 441 
(0.3614%, 1.1382%) (270, 849) 
Yukon 9.2998 % 33 
(6.5432%, 13.4933%) (23, 47) 
China 0.0002% > IM 
(0.0001%, 0.0005%) (681,427, = 1M) 
Anhui 0.0000 % Inf 
(0.0000%, 0.0000%) (Inf, Inf) 
Beijing 0.0001 % > IM 
(0.0000%, 0.0002%) (> IM,> IM) 
Chongqing 0.0000% > IM 
(0.0000%, 0.0000%) IER ERR 
Fujian 0.0000 % > IM 
(0.0000%, 0.0000%) IER a BAB) 
Gansu 0.0000 % = IM 
(0.0000%, 0.0000%) IER a BAB) 
Guangdong 0.0001 % > IM 
(0.0000%, 0.0001 %) (2 1M, = 1M) 
Guangxi 0.0005 % 676,130 
(0.0003 %, 0.0008%) (391,913, = 1M) 
Guizhou 0.0000 % > IM 
(0.0000%, 0.0000%) IER a IM) 
Hainan 0.0000 % Inf 
(0.0000%, 0.0000%) (= IM, Inf) 
Hebei 0.0000 % > IM 
(0.0000%, 0.0001%) (= IM, = IM) 
Heilongjiang 0.0004% 726,439 
(0.0002%, 0.0007%) (433,358, = 1M) 
Henan 0.0003 % 969,289 
(0.0002%, 0.0005%) (612,532, = 1M) 
Hubei 0.0000 % > IM 
(0.0000%, 0.0000%) (= 1M, = 1M) 
Hunan 0.0000 % > IM 
(0.0000%, 0.0000%) (> 1M, = 1M) 
Inner Mongolia 0.0000 % > IM 


(0.0000%, 0.0002%) 


(= 1M, = 1M) 


> IM 


(751,850, > 1M) 


939 
(578, 1,595) 
1,129 
(813, 1,919) 
531 
(360, 788) 
1,230 
(752, 2,367) 
92 
(63, 131) 
> IM 
(> IM,> IM) 
Inf 
(Inf, Inf) 
> IM 
(> IM,> IM) 
> IM 
(> IM,> IM) 
> IM 
(> IM,> IM) 
> IM 
(> IM,> IM) 
> IM 
(> IM,> IM) 
> IM 
(> IM,> IM) 
> IM 
(> IM,> IM) 
Inf 
(> IM, Inf) 
> IM 
(= 1M, = IM) 
> IM 
(> IM,> IM) 
> IM 
(> IM,> IM) 
> IM 
(> IM,> IM) 
> IM 
(= 1M, = 1M) 
> IM 
(= IM, = IM) 
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>1M 

(œ IM, = IM) 
3,596 

(2,213, 6,107) 
4,323 

(3,113, 7,346) 
2,033 

(1,379, 3,016) 
4,708 

(2,877, 9,061) 
352 

(243, 501) 

>1M 

(= 1M, = 1M) 

Inf 
(Inf, Inf) 

> IM 

(> 1M, = 1M) 
> IM 

(> IM, = IM) 
> IM 

(= 1M, = 1M) 
> IM 

(= 1M, = 1M) 
> IM 

(= 1M, = 1M) 
>1M 

(= 1M, = 1M) 
> IM 

(> IM,> IM) 

Inf 
(> IM, Inf) 

> IM 

(> IM,> IM) 
> IM 

(> IM,> IM) 
> 1M 

(> IM,> IM) 
> 1M 

(> IM,> IM) 
> IM 

(> IM,> IM) 
> IM 

(= 1M, = IM) 





> IM 

(= 1M, = 1M) 
3,801 

(2,339, 6,454) 
4,569 

(3,290, 7,764) 
2,148 

(1,458, 3,187) 
4,976 

(3,041, 9,576) 
372 

(256, 529) 

> IM 
(= 1M, = 1M) 
Inf 
(Inf, Inf) 

>1M 

(> IM, = IM) 
> 1M 

(> IM,> IM) 
> IM 

(> IM,> IM) 
> IM 

(= 1M, = 1M) 
> IM 

(> IM,> IM) 
> IM 

(= 1M, = 1M) 
> IM 

(> IM,> IM) 

Inf 
(> 1M, Inf) 

> IM 

(> IM,> IM) 
> IM 

(> IM,> IM) 
> IM 

(> IM,> IM) 
> IM 

(> IM,> IM) 
> IM 

(> IM,> IM) 
>1M 

(= 1M, = 1M) 


> IM 

(œ IM, = IM) 
3,716 

(2,287, 6,311) 
4,468 

(3,217, 7,591) 
2,101 

(1,425, 3,117) 
4,865 

(2,973, 9,363) 
364 

(251, 517) 

> 1M 

(> IM,> IM) 

Inf 
(Inf, Inf) 

> IM 

(> IM,> IM) 
> IM 

(> IM,> IM) 
> IM 

(> IM,> IM) 
> 1M 

(> IM,> IM) 
> 1M 

(> IM,> IM) 
> 1M 

(> IM,> IM) 
> 1M 

(> IM,> IM) 

Inf 
(> IM, Inf) 

> IM 

(= 1M, = 1M) 
> IM 

(= 1M, = 1M) 
> IM 

(= 1M, = 1M) 
> IM 

(> IM,> IM) 
> 1M 

(> IM,> IM) 
> IM 

(œ IM, = IM) 


> IM 

(> IM, = IM) 
2,496 

(1,536, 4,239) 
3,001 

(2,160, 5,099) 
1,411 

(957, 2,093) 
3,268 

(1,997, 6,289) 

244 
(168, 347) 
> IM 

(> IM,> IM) 

Inf 
(Inf, Inf) 

> IM 

(= 1M, = 1M) 
> 1M 

(> IM,> IM) 
> IM 

(> IM,> IM) 
> 1M 

(> IM,> IM) 
> IM 

(> IM,> IM) 
> 1M 

(> IM,> IM) 
> IM 

(> IM,> IM) 
Inf 

(> IM, Inf) 

> IM 

(> IM, = IM) 
> IM 

(= 1M, = 1M) 
> 1M 

(> IM,> IM) 
> 1M 

(> IM,> IM) 
> 1M 

(= 1M, = 1M) 
> 1M 

(œ IM, = IM) 

















Jiangsu 0.0006% 550,222 > IM > IM > IM > IM > IM 
(0.0004%, 0.0011%) (267,476, 848,138) (745,755, = 1M) (= 1M, = 1M) (> IM, = IM) (= 1M, = 1M) (= 1M, = 1M) 
Jiangxi 0.0000% > IM > IM > IM > IM > IM > IM 
(0.0000%, 0.0000%) (> IM,> IM) (> IM,> IM) (> IM,> IM) (> IM,> IM) (= 1M, = 1M) (= 1M, = 1M) 
Jilin 0.0000% > IM > IM > IM > IM > IM > IM 
(0.0000%, 0.0000%) (œ 1M, = 1M) (= 1M, = 1M) (= 1M, = 1M) (= 1M, = 1M) (= 1M, = 1M) (= 1M, = 1M) 
Liaoning 0.0016% 189,757 529,065 > IM > IM > IM > IM 
(0.0008%, 0.0028%) (111,558, 395,974) (311,038, = 1M) (= 1M, = 1M) (> IM, = IM) (= 1M, = 1M) (826,341, > 1M) 
Ningaia 0.0000% > IM > IM > IM > IM > IM > IM 
(0.0000%, 0.0001%) (> IM,> IM) (> IM,> IM) (> IM,> IM) (> IM,> IM) (> IM,> IM) (> IM,> IM) 
Qinghai 0.0000% > IM > IM > IM > IM > IM > IM 
(0.0000%, 0.0001%) (> IM,> IM) (> IM,> IM) (> IM,> IM) (> IM,> IM) (> IM,> IM) (> IM,> IM) 
Shaanxi 0.0000 % > IM > IM > IM > IM > IM > IM 
(0.0000%, 0.0001 %) (= 1M, = 1M) (> IM, = IM) (= 1M, = 1M) (= 1M, = 1M) (= 1M, = 1M) (= 1M, = 1M) 
Shandong 0.0001% > IM > IM > IM > IM > IM > IM 
(0.0001%, 0.0002%) (= 1M, = 1M) (= 1M, = 1M) (= 1M, = 1M) (> 1M, = IM) (= 1M, = 1M) (= 1M, = 1M) 
Shanghai 0.0009% 353,414 985,361 > IM > IM > IM > IM 
(0.0005%, 0.0014%) (217,455, 664,360) (606,290, > 1M) (= 1M, = 1M) (> IM, = IM) (= 1M, = 1M) (= 1M, = 1M) 
Shanxi 0.0000 % > IM > IM > IM > IM > IM > IM 
(0.0000%, 0.0000%) (> IM,> IM) (> IM, = 1M) (= 1M, = 1M) (> IM, = IM) (= 1M, = 1M) (= 1M, = 1M) 
Sichuan 0.0001% > IM > IM > IM > IM > IM > IM 
(0.0000%, 0.0002%) (= 1M, = 1M) (= 1M, = 1M) (= 1M, = 1M) (> 1M, = IM) (= 1M, = 1M) (= 1M, = 1M) 
Tianjin 0.001% 315,434 879,469 > IM > IM > IM > IM 
(0.0003%, 0.003%) (101,042, 939,120) (281,716, = 1M) (= 1M, = 1M) (= 1M, = 1M) (= 1M, = 1M) (748,440, = 1M) 
Tibet n.d. --- --- #25 === === === 
Xinjiang 0.0000% > IM > IM > IM > IM > IM > IM 
(0.0000%, 0.0000%) (> 1M, Inf) (> 1M, Inf) (> IM, Inf) (= IM, Inf) (= IM, Inf) (= IM, Inf) 
Yunnan 0.0009% 350,800 978,072 > IM > IM > IM > IM 
(0.0005%, 0.0016%) (190,158, 584,150) (530,183, = 1M) (= 1M, = 1M) (= 1M, = 1M) (= 1M, = 1M) (= 1M, = 1M) 
Zhejiang 0.0000% > IM > IM > IM > IM > IM > IM 
(0.0000%, 0.0000%) (> IM,> IM) (> IM,> IM) (> IM,> IM) (> IM, > IM) (> IM,> IM) (= 1M, = 1M) 
Egypt 0.0726% 4,226 11,781 45,099 47,661 46,602 31,300 
(0.0273 %, 0.1138%) (2,697, 11,256) (7,518, 31,382) (28,780, 120,129) (30,415, 126,955) (29,740, 124,133) (19,974, 83,373) 
European Union 3.9648 % 71 216 826 873 854 573 
(3.0069%, 5.4727%) (56, 102) (156, 285) (598, 1,089) (632, 1,151) (618, 1,125) (415, 756) 
Austria 7.582% 40 113 432 456 446 300 
(7.2889%, 9.6649%) (32, 42) (89, 117) (339, 449) (358, 475) (350, 464) (235, 312) 
Belgium 6.6396% 46 129 493 521 510 342 
(4.5097 %, 8.7201%) (35, 68) (98, 190) (376, 726) (397, 767) (388, 750) (261, 504) 
Bulgaria 2.2831% 134 375 1,434 1,516 1,482 996 
(1.5615%, 3.1896%) (96, 197) (268, 548) (1,027, 2,097) (1,085, 2,217) (1,061, 2,167) (713, 1,456) 
Croatia 6.2509% 49 137 524 554 541 364 
(5.4277%, 7.3272%) (42, 57) (117, 158) (447, 603) (472, 638) (462, 623) (310, 419) 
Cyprus 1.8986% 162 451 1725 1,823 1,782 1,197 
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(1.3388%, 2.5765%) (119, 229) (332, 639) (882, 1,698) 
Czech Republic 10.128% 30 84 323 342 334 224 
(9.0398%, 10.9588%) (28, 34) (78, 95) (299, 362) (316, 383) (309, 374) (207, 251) 
Denmark 4.4201% 69 194 741 783 766 514 
(3.237%, 5.9783%) (51, 95) (143, 264) (548, 1,012) (579, 1,069) (566, 1,045) (380, 702) 
Estonia 3.2309% 95 265 1,014 1,071 1,047 703 
(2.4958%, 4.555%) (67, 123) (188, 343) (719, 1,312) (760, 1,387) (743, 1,356) (499, 911) 
Finland 1.7646% 174 485 1,856 1,961 1,918 1,288 
(1.0067%, 3.4921%) (88, 305) (245, 850) (938, 3,253) (991, 3,438) (969, 3,361) (651, 2,258) 
France 0.8692% 353 984 
(0.3183%, 1.6579%) (185, 964) (516, 2,688) 
Germany 2.9448 % 104 291 1,112 1,175 1,149 772 
(2.2657%, 7.0048%) (44, 135) (122, 378) (468, 1,445) (494, 1,528) (483, 1,494) (324, 1,003) 
Greece 3.4772% 88 246 942 995 973 654 
(2.5819%, 5.7358%) (53, 119) (149, 331) (571, 1,268) (603, 1,340) (590, 1,311) (396, 880) 
Hungary 5.4027% 57 158 606 641 626 421 
(3.4865 %, 8.0329%) (38, 88) (107, 245) (408, 939) (431, 993) (421, 971) (283, 652) 
Ireland 2.7065% 113 316 1.210 1,279 1,250 840 
(1.3036%, 4.41%) (70, 235) (194, 656) (743, 2,512) (785, 2,655) (767, 2,596) (515, 1,744) 
Italy 0.8517% 360 1,005 
(0.6251%, 1.2349%) (248, 491) (693, 1,369) 
Latvia 4.7632% 64 180 688 727 710 477 
(3.2874%, 6.8327%) (45, 93) (125, 260) (479, 996) (507, 1,053) (495, 1,029) (333, 691) 
Lithuania 3.1464% 98 272 1,041 1,100 1,076 722 
(2.4847%, 4.4227%) (69, 123) (193, 344) (740, 1,318) (783, 1,393) (765, 1,362) (514, 915) 
Luxembourg 1.7269% 178 495 1,316 
(1.211%, 2.6969%) (114, 253) (317, 706) (843, 1,877) 
Malta 0.9345 % 328 915 
(0.6545%, 1.2744%) (241, 469) (671, 1,307) 
Netherlands 6.4673 % 47 132 506 523 351 
(4.752%, 8.3297%) (37, 65) (103, 180) (393, 689) (406, 712) (273, 478) 
Poland 2.118% 145 404 1,073 
(1.9921%, 2.6347%) (116, 154) (325, 429) (863, 1,141) 
Portugal 0.5193% 591 
(0.2801%, 0.7387%) (415, 1,096) 
Romania 4.4927% 68 190 729 770 753 506 
(2.8969%, 6.0855%) (50, 106) (141, 295) (538, 1,131) (569, 1,195) (556, 1,168) (373, 785) 
Slovakia 9.1138% 34 94 359 380 371 249 
(7.8365%, 10.6227%) (29, 39) (81, 109) (308, 418) (326, 442) (319, 432) (214, 290) 
Slovenia 9.2961% 33 92 352 372 364 244 
(6.2612%, 14.0264%) (22, 49) (61, 137) (233, 523) (247, 553) (241, 540) (162, 363) 
Spain 0.0575% 
(0.0369%, 0.0874%) 
Sweden n.d. --- --- --- --- --- --- 
Iceland 2.197% 140 389 1,491 1,575 1,540 1,035 
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(1.3208%, 4.0252%) (76, 232) (213, 648) (814, 2,480) (860, 2,620) (841, 2,562) (565, 1,721) 
India 0.0904% 3,395 9,466 36,235 38,294 37,443 25,148 
(0.0789%, 0.2077%) (1,477, 3,887) (4,119, 10,836) (15,769, 41,482) (16,665, 43,838) (16,294, 42,864) (10,944, 28,789) 
Indonesia 0.0101% 30,352 84,624 
(0.005%, 0.0217%) (14,139, 61,685) (39,422, 171,986) 
Israel 0.1221% 2,513 7,007 26,821 28,345 27,715 18,615 
(0.0792%, 0.1805%) (1,700, 3,876) (4,740, 10,806) (18,144, 41,365) (19,175, 43,716) (18,749, 42,744) (12,592, 28,709) 
Japan 0.0056% 55,018 153,397 
(0.0034%, 0.0119%) (25,686, 90,788) (71,615, 253,127) 
Mexico 0.2971 % 1,033 2,879 11,022 11,648 11,389 7,649 
(0.1372%, 0.68%) (451, 2,236) (1,258, 6,235) (4,816, 23,867) (5,090, 25,223) (4,977, 24,663) (3,343, 16,565) 
New Zealand 0.1437% 2,136 5,955 22,795 24,090 23,555 15,820 
(0.1056%, 0.2102%) (1,460, 2,907) (4,070, 8,105) (15,578, 31,027 (16,463, 32,790 (16,097, 32,061) 10,812, 21,534) 
Nigeria 0.0051% 60,720 169,295 
(0.0025%, 0.0109%) (28,140, 124,082) (78,459, 345,955) 
Norway 2.7376% 112 313 1,196 1,264 1,236 830 
(2.0521%, 4.0047%) (77, 150) (214, 417) (818, 1,596) (864, 1,687) (845, 1,649) (568, 1,108) 
Qatar 0.4208 % 729 2,033 7,783 8,226 8,043 5,402 
(0.236%, 0.5729%) (536, 1,300) (1,493, 3,626) (5,716, 13,879) (6,041, 14,668) (5,907, 14,342) (3,967, 9,633) 
Russia 0.7135% 430 1,199 4,590 4,850 4,743 3,185 
(0.6167%, 0.8016%) (383, 498) (1,067, 1,387) (4,086, 5,310) (4,318, 5,612) (4,222, 5,487) (2,836, 3,686) 
Saudi Arabia 0.002% 150,840 
(0.0006%, 0.0052%) (59,130, 482,837) 
Serbia 1.3359% 230 640 2,451 2,591 2,533 1,701 
(0.9266%, 2.0453%) (150, 331) (418, 923) (1,601, 3,534) (1,692, 3,735) (1,655, 3,652) (1,111, 2,453) 
South Africa 0.048% 6,399 17,841 68,293 72,173 70,569 47,397 
(0.025%, 0.1222%) (2,511, 12,254) (7,002, 34,164) (26,803, 130,780) (28,326, 138,210) (27,697, 135,139) (18,602, 90,765) 
South Korea 0.3687% 832 2,320) 8,881 9,386 9,177 6,164 
(0.2869%, 0.6332%) (485, 1,069) (1,351, 2,982) (5,172, 11,414) (5,466, 12,063) (5,345, 11,795) (3,590, 7,922) 
Switzerland 1.2712% 241 673 2,576 2,723 2,662 1,788 
(0.9103%, 2.7461%) (112, 337) (312, 940) (1,193, 3,597) (1,260, 3,802) (1,232, 3,717) (828, 2,497) 
Turkey 0.5889% 521 1,453 5,561 5,877 5,747 3,860 
(0.495%, 0.6974%) (440, 620) (1,227, 1,728) (4,696, 6,616) (4,962, 6,992) (4,852, 6,836) (3,259, 4,592) 
Ukraine 4.3754% 70 196 748 791 773 519 
(2.9122%, 5.3887%) (57, 105) (159, 294) (608, 1,125) (642, 1,188) (628, 1,162) (422, 780) 
United Arab 0.094% 3,265 9,103 34,847 36,827 36,009 24,185 
Emirates (0.0494%, 0.2498%) (1,229, 6,211) (3,426, 17,316) (13,113, 66,286) (13,858, 70,053) (13,550, 68,496) (9,101, 46,005) 
United Kingdom 2.5715% 119 333 1,274 1,346 1,316 884 
(2.0297%, 3.4445%) (89, 151) (248, 421) (951, 1,613) (1,005, 1,705) (982, 1,667) (660, 1,120) 
United States of 1.7919% 171 477 1,828 1,931 1,889 1,268 
America (1.1619%, 2.687%) (114, 264) (318, 736) (1,219, 2,819) (1,288, 2,979) (1,259, 2,912) (846, 1,956) 
Alabama 1.6329% 188 524 2,006 2,120 2,072 1,392 
(0.9706%, 2.5944%) (118, 316) (330, 881) (1,262, 3,374) (1,334, 3,566) (1,304, 3,487) (876, 2,342) 
Alaska 3.9535% 78 216 828 875 856 575 
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Arizona 
Arkansas 
California 
Colorado 
Connecticut 
Delaware 
District of Columbia 
Florida 
Georgia 
Hawaii 
Idaho 
Illinois 
Indiana 
Iowa 
Kansas 
Kentucky 
Louisiana 
Maine 
Maryland 
Massachusetts 
Michigan 


Minnesota 


(2.5054%, 5.87%) 
2.7643% 
(1.8833%, 3.7897%) 
0.5864% 
(0.344%, 0.9375%) 
1.9291% 
(1.1926%, 2.9288%) 
3.8648% 
(2.5974%, 5.5929%) 
0.7765% 
(0.3348%, 1.3932%) 
1.4611% 
(0.8877%, 2.1985%) 
0.9928% 
(0.4965%, 1.8429%) 
0.1491% 
(0.0735%, 0.3026%) 
0.3715% 
(0.1839%, 0.7033%) 
1.1000% 
(0.6391%, 1.6064%) 
1.7874% 
(1.1244%, 2.779%) 
1.6601% 
(1.034%, 2.4965%) 
1.887% 
(1.2236%, 2.7637%) 
3.1211% 
(2.073%, 4.6638%) 
2.205% 
(1.4071%, 3.3584%) 
1.3104% 
(0.7906%, 2.1014%) 
0.1816% 
(0.0994%, 0.3329%) 
4.8815% 
(3.4574%, 6.8067%) 
1.2408% 
(0.7756%, 1.9289%) 
1.5626% 
(0.995%, 2.3719%) 
3.122% 
(2.3887%, 4.1254%) 
3.7262% 
(2.5439%, 5.2348%) 





(52, 122) 
111 
(81, 163) 
523 
(327, 892) 
159 
(105, 257) 
79 
(55, 118) 
395 
(220, 917) 
210 
(140, 346) 
309 
(166, 618) 


826 
(436, 1,668) 
279 
(191, 480) 
172 
(110, 273) 
185 
(123, 297) 
163 
(111, 251) 


(91, 218) 
234 
(146, 388) 
1,689 
(922, 3,086) 


(159, 396) 
196 
(129, 308) 
98 
(74, 128) 
82 
(59, 121) 



























(146, 341) (558, 1,307) (590, 1,381) (577, 1,351) (387, 907) 
309 1,185 1,252 1,224 822 
(226, 454) (864, 1,739) (913, 1,838) (893, 1,797) (600, 1,207) 
1,459 
(913, 2,487) 

443 1,178 
(292, 717) (776, 1,906) 
221 847 896 876 588 
(153, 329) (586, 1,261) (619, 1,333) (605, 1,303) (406, 875) 

1,102 
(614, 2,555) 
586 
(389, 964) 
862 
(464, 1,723) 


718 














(533, 1,339) 

479 1,272 
(308, 761) (818, 2,021) 
515 1,369 
(343, 827) (910, 2,198) 
453 1,204 
(310, 699) (822, 1,858) 
274 1,049 1,109 1,084 728 
(183, 413) (702, 1,580) (742, 1,670) (726, 1,633) (487, 1,096) 
388 1,485 1,031 
(255, 608) (975, 2,327) (677, 1,615) 


653 








(407, 1,082) 
175 671 709 693 466 
(126, 247) (481, 947) (508, 1,001) (497, 979) (334, 657) 
689 
(444, 1,103) 
547 1,455 
(361, 860) (958, 2,284) 
274 1,049 1,109 1,084 728 
(207, 358) (794, 1,371) (839, 1,449) (820, 1,417) (551, 952) 
230 879 929 908 610 
(163, 336) (626, 1,287) (661, 1,360) (646, 1,330) (434, 893) 
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Mississippi 
Missouri 
Montana 
Nebraska 
Nevada 

New Hampshire 
New Jersey 
New Mexico 
New York 
North Carolina 
North Dakota 
Ohio 
Oklahoma 
Oregon 
Pennsylvania 
Rhode Island 
South Carolina 
South Dakota 
Tennessee 
Texas 

Utah 

Vermont 


Virginia 


0.3868% 
(0.2061%, 0.6979%) 
1.9784% 
(1.3369%, 2.8753%) 
2.493% 
(1.5002%, 3.9065%) 
2.8688% 
(1.7603%, 4.6931%) 
1.6236% 
(0.9193%, 2.587%) 
5.5646% 
(3.4291%, 7.8332%) 
1.3089% 
(0.9035%, 1.805%) 
4.3246% 
(2.8308%, 6.1917%) 
3.2186% 
(2.3387%, 4.3285%) 
0.8165% 
(0.4574%, 1.4275%) 
2.5186% 
(1.6199%, 3.8533%) 
2.4753% 
(1.7355%, 3.4126%) 
1.11% 
(0.6703%, 1.785%) 
1.938% 
(1.114%, 3.1602%) 
3.1145% 
(2.2874%, 4.1416%) 
1.3159% 
(0.8309%, 1.9954%) 
0.3772% 
(0.1976%, 0.7019%) 
1.9808% 
(1.2039%, 3.1375%) 
0.5464% 
(0.2905%, 0.9501 %) 
0.4587% 
(0.2479%, 0.8396%) 
1.3594% 
(0.906%, 1.9593%) 
5.7049% 
(4.053%, 7.6921%) 
1.5061% 


793 
(440, 1,489) 
155 
(107, 230) 
123 
(19, 205) 
107 
(65, 174) 
189 
(119, 334) 
55 
(39, 89) 
234 
(170, 340) 
71 
(50, 108) 


(215, 671) 
122 
(80, 189) 
124 
(90, 177) 
276 
(172, 458) 
158 
(97, 275) 
99 
(14, 134) 
233 
(154, 369) 
814 
(437, 1,553) 
155 
(98, 255) 
562 
(323, 1,056) 
669 
(365, 1,238) 
226 
(157, 339) 
54 
(40, 76) 
204 










432 1,149 

































(298, 640) (790, 1,700) 
343 1,314 1,388 1,357 912 
(219, 570) (838, 2,183) (886, 2,307) (866, 2,256) (582, 1,515) 
298 1,142 1,206 1,180 792 
(182, 486) (698, 1,860) (737, 1,966) (721, 1,922) (484, 1,291) 
527 1,400 
(331, 931) (879, 2,472) 
154 589 622 608 408 
(109, 249) (418, 955) (442, 1,009) (432, 987) (290, 663) 
654 
(474, 947) 

198 757 800 783 526 
(138, 302) (529, 1,157) (559, 1,223) (547, 1,195) (367, 803) 
266 1,018 1,075 1,051 706 
(198, 366) (757, 1,400) (800, 1,480) (782, 1,447) (525, 972) 
1,048 
(599, 1,870) 

340 1,300 1,374 1,344 902 
(222, 528) (850, 2,022) (898, 2,137) (878, 2,089) (590, 1,403) 
346 1,323 1,367 918 
(251, 493) (960, 1,887) (992, 1,950) (666, 1,310) 
771 
(479, 1,276) 

441 1,173 
(271, 768) (719, 2,040) 
275 1,052 1,111 1,087 730 
(207, 374) (791, 1,432) (836, 1,513) (817, 1,479) (549, 994) 
650 
(429, 1,030) 











432 
(273, 711) 
1,566 
(900, 2,945) 


1,147 
(724, 1,888) 





629 


(437, 944) 
150 574 607 593 398 
(111, 211) (426, 808) (450, 854) (440, 835) (295, 561) 
568 1,509 
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(0.8939%, 2.369%) (130, 343) (361, 957) (959, 2,543) 
Washington 1.5935% 193 537 1,426 
(0.8836%, 2.6517%) (116, 347) (323, 968) (857, 2,572) 
West Virginia 2.2764% 135 376 1,439 998 
(1.5275%, 3.2917%) (93, 201) (260, 560) (995, 2,144) (690, 1,488) 
Wisconsin 3.3275% 92 257 984 1,040 1,017 683 
(2.256%, 4.7204%) (65, 136) (181, 379) (694, 1,452) (733, 1,534) (717, 1,500) (482, 1,007) 
Wyoming 2.2208% 138 385 1,475 1,023 
(1.4641%, 3.2833%) (93, 210) (261, 584) (997, 2,237) (692, 1,552) 











95% CI=95% confidence interval. Brackets are the 95% confidence intervals (lower limit, upper limit). n.d.=no data. NNE=number needed to exclude. NNEs are color-coded by the lower 
limit of the 95% CI: yellow (= 1,000 and < 10,000), orange (= 10,000 and < 100,000), red (= 100,000 and < 1,00,000), and black (= 1,000,000). NNEs = 1,000,000 were truncated to > 1M. 
+Estimates of the baseline infection risk (i.e., point-prevalence of infectious cases in the jurisdiction) were taken from a database (https://decision-support-tools.com/) developed and regularly 
updated by Defence Research and Development Canada (DRDC), an agency of the Canadian Department of National Defence. Baseline infection risks are shown at most to the fourth decimal 
place due to space limitations. A true 0% baseline infection risk can be identified in the table when the NNE is infinite (Inf) because that means there were no infectious cases 
(NNE=1/0=infinite). 
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DISCUSSION 

There were two main findings of our study. First, the NNEs suggest that at least 1,000 
unvaccinated people likely need to be excluded to prevent one SARS-CoV-2 transmission event 
in most types of settings for many jurisdictions, notably Australia, California, Canada, China, 
France, Israel, and others. Second, the NNEs of almost every jurisdiction examined are well 
within the range of the NNTs of ASA in primary prevention of CVD (= 250 to 333). There are 
four main implications of our study. 

First, the two conservative assumptions of the model mean that the calculated NNEs 
likely represent the minimum number of unvaccinated people needed to exclude to prevent one 
SARS-CoV-2 transmission event in each jurisdiction and type of setting. The true NNEs are 
likely much higher, especially in jurisdictions with higher levels of vaccine/natural immunity. 

Second, while SARS-CoV-2 vaccines are beneficial, the high NNEs suggest that 
excluding unvaccinated people has negligible benefits for reducing SARS-CoV-2 transmissions 
in many jurisdictions across the globe. This is because unvaccinated people are likely not at 
significant risk — in absolute terms — of transmitting SARS-CoV-2 to others in most types of 
settings (as of mid-to-end November 2021). This is why so many unvaccinated people likely 
need to be excluded to prevent one transmission event. While a meta-analysis found that the 
RRR of vaccines on SARs is ~41%,° current baseline transmission risks in most settings are 
negligible, such that a 41% decrease of a negligible risk is still a negligible absolute risk 
difference between vaccinated vs. unvaccinated people. This does not mean vaccines are not 
beneficial for reducing transmission; on the contrary, a 41% RRR in SARs in a population with 
high vaccination levels will lower the burden of SARS-CoV-2. Rather, it means that excluding 


unvaccinated people on the basis of the relative risk difference between vaccinated vs. 
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unvaccinated people is likely not justifiable because baseline transmission risks are negligible for 
most types of settings and jurisdictions. 

Third, societies and policymakers should urgently consider the harms of exclusion via 
VMVP. This is because these NNEs suggest that exclusion may not be a proportionate response 
to the risk unvaccinated people pose to others. Just as ASA is not recommended for primary 
prevention of CVD because the harms outweigh the benefits (NNTs = 250 to 333), the harms of 
excluding unvaccinated people via VMVP may outweigh the benefits (NNEs > 250 to 333). The 
analogy to ASA is made not because the outcomes are the same. Clearly, they are not: one is a 
within-individual risk (e.g., myocardial infarction), the other is a between-individual or 
‘externalized’ risk where one or more other people may be impacted (i.e., transmission event). 
Rather, we are drawing a comparison between a public health intervention (i.e., exclusion) and 
an intervention widely recognized in medicine as having limited benefit due to the negligible 
ARRs. Thus, while the outcomes are different, the NNTs of ASA in primary prevention of CVD 
provide a benchmark for interpreting the magnitude of the NNEs. Based on this benchmark, our 
study suggests that this public health intervention (exclusion) has limited benefit. Given that 
many jurisdictions have or are planning to implement exclusionary policies via VMVP, there is 
an urgent need to consider the harms in light of the limited benefits. At least four harms should 
be considered. First, there is the possibility of staffing shortages from the loss of unvaccinated 
healthcare workers, which may have downstream effects on patient care. Second, there is the 
harm of unemployment and — in some cases — unemployability if vaccines are a condition of 
employment or higher-education in most organizations. Third, this would mean unvaccinated 


people and their families will face financial hardship. This may have downstream effects on their 
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mental/physical health. A fourth harm of exclusion is the creation of a class of citizens who are 
not allowed to fully participate in many areas of society. 

Fourth, we suggest that the NNE and the baseline infection risk Oe, the point-prevalence 
of infectious cases) be included as key indicators for determining the harms vs. benefits of 
VMVP. It compliments other indicators, such as hospital/ICU occupancy, the basic and effective 
reproduction number, SARS-CoV-2-related hospitalizations/deaths, and regional vaccination 
rates. Moreover, because the model can be updated to account for any future potentially 
dominant variants (e.g., Omicron), the NNE enables efficient monitoring of the benefits of 
exclusion in real-time as the pandemic changes and variants emerge. This will improve public 
health because the NNE provides societies and policymakers with an intuitive metric to monitor 
whether to implement exclusionary policies via VMVP and when to stop these policies. 

These findings depend on the accuracy of the assumptions of the NNE model and the 
accuracy of baseline infection risks estimates. The first assumption was that the Delta variant 
correction factor (1.97) we applied to the robust estimates of the wild-type SARs was accurate. 
We discovered that this was an accurate but conservative assumption since the corrected robust 
SAR in households was 26% higher than what has been observed in real-world transmission 
studies of the Delta variant in unvaccinated index cases.”**° The higher corrected SARs 
(SARw, X CFpeita) of the model increased the estimated baseline transmission risks, which 
lowered the NNEs. The second assumption was that the SARs for each type of setting are 
relatively consistent. This was supported by using robust estimates of the SARs from multiple 
meta-analyses. The third assumption was that the baseline infection risk is relatively stable 
across individuals in a population. This was addressed by calculating 95% Cls around the NNEs 


for each jurisdiction using the 95% CIs around the baseline infection risks. The fourth 
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assumption was that (i) unvaccinated people have no natural immunity and (ii) the contacts of 
infected unvaccinated index cases have no vaccine or natural immunity. This was a strong 
conservative assumption because it means the estimated baseline transmission risks were higher 
than what would be predicted in a population with some level of vaccine/natural immunity 
because studies show that vaccine/natural immunity reduces SARs.” 23 

Estimating the baseline infection risk in a population is challenging because it requires 
modeling various viral factors (e.g., incubation and infectious periods, asymptomatic cases) in 
addition to case counts to arrive at the current point-prevalence of infectious cases. There is 
likely some degree of underestimation of the baseline risk in some jurisdictions due to 
underreporting of cases. However, it is extremely difficult to accurately model the 
underreporting rate for each jurisdiction and how it changes over time because one is trying to 
model something where there is no data. This is shown by the extreme variation in 
underreporting estimates.”!?” Local context/knowledge remains essential for estimating 
underreporting rates in a jurisdiction over time. Since this is not available at a global scale, the 
degree of underestimation is unknown. Underestimation was partially addressed by utilizing the 
95% CIs of the baseline infection risks in order to represent the likely general range of the 
baseline infection risks. The conservative assumptions also partially compensate for this because, 
as we just noted, they increased the estimated baseline transmission risks, which is what 
correcting for underreporting would do. Notwithstanding this limitation, we have provided a 
‘proof of concept’ of how baseline infection risk data can be used to calculate the NNE to 
evaluate the benefits of excluding unvaccinated people to reduce SARS-CoV-2 transmissions. 

Future research should focus on developing an online interactive global database of the 


NNEs and baseline infection risks in each region which is regularly updated. Furthermore, this 
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interactive database should be linked to regional vaccination rates, SARS-CoV-2-related 
hospitalizations/deaths, hospital/ICU occupancy, and the basic and effective reproduction 
number to contextualize the NNEs within other important metrics. This would be helpful for 
societies and policymakers to evaluate in real-time the harms vs. benefits of excluding 
unvaccinated people. Health economic analyses are also needed to analyze the cost-effectiveness 
and cost-benefits of excluding unvaccinated people. Future studies could incorporate 
vaccine/natural immunity data into the model to explore the impact of herd immunity on the 
NNE. Finally, when new variants emerge (e.g., Omicron) and potentially outcompete the Delta 


variant, the NNEs should be updated with the transmissibility data of future dominant variants. 
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